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ABSTRACT 

We study a 24 /mi selected sample of 330 galaxies observed with the Infrared 
Spectrograph for the 5 mJy Unbiased Spitzer Extragalactic Survey. We estimate 
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accurate total infrared luminosities by combining mid-IR spectroscopy and mid- 
to-far infrared photometry, and by utilizing new empirical spectral templates 
from Spitzer data. The infrared luminosities of this sample range mostly from 
1O 9 L to 1O 13 - 5 L , with 83% in the range 1O 1O L <L IR <1O 12 L . The redshifts 
range from 0.008 to 4.27, with a median of 0.144. The equivalent widths of the 
6.2 fim aromatic feature have a bimodal distribution. We use the 6.2 /im PAH 
EW to classify our objects as SB-dominated (44%), SB-AGN composite (22%), 
and AGN-dominated (34%). The high EW objects (SB-dominated) tend to have 
steeper mid-IR to far-IR spectral slopes and lower Lfr and redshifts. The low 
EW objects (AGN-dominated) tend to have less steep spectral slopes and higher 
Lir and redshifts. This dichotomy leads to a gross correlation between EW and 
slope, which does not hold within either group. AGN dominated sources tend 
to have lower log(LpAH7.7/jm/LpAHii.3/jm) ratios than star-forming galaxies, pos- 
sibly due to preferential destruction of the smaller aromatics by the AGN. The 
log(LpAH7.7^m/LpAHii.3^m) ratios for star-forming galaxies are lower in our sam- 
ple than the ratios measured from the nuclear spectra of nearby normal galaxies, 
most probably indicating a difference in the ionization state or grain size distri- 
bution between the nuclear regions and the entire galaxy. Finally, we provide 
a calibration relating the monochromatic 5.8, 8, 14 and 24 /mi continuum or 
Aromatic Feature luminosity to Lip for different types of objects. 

Subject headings: galaxies: active, galaxies: starburst, galaxies: evolution, in- 
frared radiation, surveys 



1. Introduction 



Infrared bright galaxies play critical roles in galaxy formation and evolution. The In- 
fraRed Astronomical Satellite (IRAS) facilitated the s tudy of an import a nt group of objects, 
the U ltra Luminous InfraRed Galaxies (ULIRGs) (jSoifer et al.l Il989t ISanders fc Mirabel 



19961 ) , which were first hinted at by ground base d observations of iRieke fc Low! ( 119721 ) . 
Studies fro m the Infrared Spa c e Observatory (I SO) (jElbaz et al.lll999l ) and the Spitzer Space 
Telescope (IHouck et al.l 120051 : lYan et al.l 120071 ) later revealed that LIRGs and ULIRGs are 
much more common at high redshift than in the local Universe. The number density o f 
IR luminous galaxies evolves strongly with redshift to at least z~l (|Le Floc'h et al.l 120051 ). 
The fraction of galaxies powered by star formation versus AGN is still controversial, but is 
crucial for determining unbiased luminosity functions for various categories of objects and 
understanding the evolution process. 
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The superb sensitivity of the Spitzer Space Telescope ( IWerner et al.ll2004l ) has led to the 
discovery of new populations of faint, high-redshift galaxies with extreme IR/optical colors 
jDickinson et al.lboO^ Irlouck et al.lbood: IWeedman et al-l iooS lYan et al.l2007t ICaputi et al~ 



20071 : iDev et al 



2008 



Dasyra et al.l 120091 ). However, these studies often have at least 



one other constraint than the mid-IR flux limit, usually a minimum R band magnitude 
or an IRAC-based color selection, designed to favor sources in specific redshift ranges, or 
with high luminosity. The 5 Millijanksy Unbiased Spitzer Extragalactic Survey (5MUSES) 
is an infrared selected sample. A major advantage of 5MUSES is its simple selection: 
f t/ (24/im)>5mJy. This relatively bright flux limit allows for a more detailed study of the 
infrared properties, filling in the gap between local galaxies and high redshift samples, and 
helping to improve the modeling of galaxy populations and their evolution. 

In order to advance our understanding of the properties and evolution of galaxies, it 
is crucial to obtain accurate estimates of their bolometric luminosities. Several studies 
have shown that monochromatic lumino s ities in the mid-IR can be used to est imate Lpr 
JSaiina et al.l I2OO7I : iBavouzet etaD 120081 : iRieke et al.l I2OO9I : ICalzetti et all boioh . and the 
uncertaint ies on these estimates decrease significantly when far-infrared (FIR) fluxes are 
available (IKartaltepe et al.ll2010l ). However, the spectral energy distr ibution (SEP) of star- 



forming galaxies, AGN and ULIRGs display a wide range of shapes ( Weedman et al. 



Brandl et al. 2006: Smith et al. 2007: Armus et al. 2007: Haoet al. 2007; Wu et al 



2005; 



20091 : 

Veilleux et al.l 120091 ) . Applying these methods without knowing a source's spectral type 
could cause significant biases in luminosity estimates between types of objects and seriously 
mislea d the interpretatio ns. The 5 to 36 fim spectra obtained by the Infrared Spectrograph 
(IRS) (IHouck et al.ll2004j ) for the 5MUSES sample allows for aromatic feature identification, 
excitation line analysis, and decomposition into star formation and AGN components, thus 
providing essential information for classifying the origin of the luminosity. 

The mid-IR is home to a set of broad emissio n line features, which are thought to 



origin ate from Polycyclic Aromatic Hydrocarbons (IPuget et al.l Il985l ; lAUamandola et al. 
1989). PAHs are organic molecules that are ubiquitous in o ur own Galaxy (jPeeters et al.l 



20021 ) and nearby star-forming galaxies (IHelou et al.l l200ll ; ISmith et al.l 120071 ) . In total 
they can contribute a significant fraction (10% or more) of the tota l infrared luminosity 
in star-forming 



Wu et al. 



2006; 



galaxies. PAHs are wea k in low metallicity galaxies dMadden et al. 



Weedman et al 



Engelbracht et al 



2005 



AGN ( Smith et al. 



2008|), or in galaxies wi t h powerful (|Roche et al. 



Armus et al. 



2007 



Dale et al. 



20071 : iDesai et all 120071 : IWu et al.l 120091 ) or even weak 



2006 



1991 



20091 ). The PAH features, including their p rofiles, central 



wavelengths and band-to-b and intensity ratios have been s tudied i n deta il by iPeeters et al. 



(120021 ). ISmith et al.l (120071 ) and most recently reviewed by iTielend (120081 ) . The 6.2 /im fea- 
ture and the 7.7 /xm complex are attributed to vibrational modes of the carbon skeleton. 
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The 8.6 /zm feature is attributed to in-plane C-H bending, while the features at 11.3/im 
and 12.7 /zm are identified as out-of-plane C-H bending modes. It is generally thought that 
charged PAHs radiate more strongly in the C-C vibrational modes, while neutral PAHs radi- 
ate strongly in the out-of-plane C-H bending modes at 11.3 fim and 12.7/xm. The fraction of 
the power radiated by PAH in the different bands followi ng single-photon h eating depends 
on both the PAH ionization and on the size of the PAH (jDraine fc Lil 120071 ). Thus the ob- 
served variations in the PAH band-to-band ratios can reflect yariations in physica l conditions 
JSmith et al. l l2007l : iGalliano et alll2008[ bordon et aDbood : b'Dowd et aDl20Q9h . 



Because PAH emission can be very prominent in star-forming systems, it has often been 
used as a relatively extinction-free diagnostic to ol to constrain star formation. Detailed 



studies on the pro perties of PAH features locally (ISpoon et al. 



2007 



Desai et al.l 120071 ) and 



bpc 

at higher redshift (lYan et al.ll2005l ; iHouck et al.ll2005l ; iHuang et al.ll2009l ) reveal differences in 
the PAH Equivalength Widths (EWs) and Lpah/Lir ratios. This might indicate that some 
evolution in the PAH properties occurs with redshift, or that sample selection effects make 
for large variations in the aromatic feature properties. However, one cannot simply apply 
our knowledge from the local universe to high redshift galaxies, or make fair comparisons 
between the two unless truly equivalent samples have been studied. Current analysis on the 
PAH properties are based on ISO or Spitzer observations of relatively bright objects, which 
have been selected because of previously known optical or IRAS criteria. Thus it is crucial 
to have a complete or at least unbiased census of galaxies in order to understand the galaxy 
evolution process and its relation to the aromatic feature emission. 

In this paper, we study the properties of PAH emission and IR luminosities. This is the 
first of a series of papers to study the IR selected representative sample of 5MUSES. Helou 
et al. (in preparation) will address the general properties of the sample and how it bridges 
the gap between local and high-z galaxies. Yong et al. (in preparation) will present the 
correlations between old stars and current star formation. Detailed population modelling 
will also be performed to address the bimodal distribution of the PAH EWs discovered 
in this study. In Section 2, we briefly describe the sample selection, data reduction and 
measurements of spectral features. We introduce our library of empirical IR SED templates 
built upon Spitzer observations in Section 3, and derive the total infrared luminosities for 
5MUSES galaxies. We also discuss how well one can constrain the IR SED if only mid-IR 
data are available. In Section 4, we study the properties of PAH emission from our flux 
limited sample. Finally, we present our conclusions in Section 5. Using the IR luminosities 
we derived in Section 3 and the PAH luminosities from Section 4, we discuss estimation of Lir 
from PAH luminosity or monochromatic continuum luminosity in the Appendix. Throughout 
this work, we assume a ACDM cosmology with Ho=70kms~ 1 Mpc -1 , f2 m =0.27 and Q\=0.73. 
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2. Observations and Data Analysis 
2.1. The Sample 



5MUSES is a mid-IR spectroscopic survey of a 24 /im flux-limited (f24^m >5mJy) 
representative sample of 330 galaxies. The galaxies are selected from the SWIRE fields 
(Lonsdale et~aD l2003h . including Elais-Nl (9.5 deg 2 ), Elais-N2 (5.3 deg 2 ), Lockman Hole 
(11.6 deg 2 ) and XMM (9.2 deg 2 ), in addition to the Spitzer Extragalactic First Look Sur- 
vey (XFLS, 5.0 deg 2 ) field (IFadda et al.ll2006l ). It provides a representative sample at in- 
termediate redshi: 
forming galaxies ( 



t (< z >~0.144) which bridges the gap between the brig ht, nearby star 



Kennicutt et al.1 2003 



Smith et all 120071 bale et all 120091 ). local ULIRGs 



( lArmus et al.ll2007l ; iDesai et al.ll2007l ; IVeilleux et al .11200 91 ) and t he much fainter and more dis 
tant sources pursued in most z~2 IRS follow-up work to date (IHouck et al.ll2005l ; lYan et al. 
20071 ). The full details of the sample, including selection criteria and observation strategy 
are covered in Helou et al. (2010, in prep). 



2.2. Observation and Data Reduction 

Because of its selection in the SWIRE and XFLS fields, IRAC 3.6-8.0 //m photometry 
is available for the entire 5MUSES sample. In addition to the MIPS 24 /im photometry used 
to select this sample, 90% of our sources have also been detected at MIPS 70 /im and 54% 
have been detected at MIPS 160 /im. Low-resolution spectra (R=64~128) of all 330 galaxies 
in 5MUSES have been obtained with the Short-Low (SL: 5.5-14.5 /im and Long-Low (LL: 
14-35 /im) modules of the IRS using the staring mode observations. The integration time 
on each object was estimated based on its 24 /im flux densities and typically ranges from 
300-960 seconds (see Tabled]). A sub-set of the 5MUSES sample has also been observed 
with the high-resolution modules of the IRS, which will be covered in a future paper. 

The low-resolution IRS data were processed by the Spitzer Science Center data reduc- 
tion pipeline version SI 7. The two-dimensional image data were converted to slopes after 
linearization correction, subtraction of darks, cosmic-ray removal, stray light and flat field 
correction. The post-pipeline reduction of the spectral data started from the pipeline prod- 
ucts basic calibrated data (bed) files. We took the median of all images from the off-source 
part of the slit (off-order and off-nod) and then subtracted it from the image on the source. 
Then we combined all the background-subtracted images at one nod position and took the 
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mean. The resulting images were then cleaned with the IRSCLEAN package^ to remove bad 
pixels and apply rogue pixel correction. 

We used the Spitzer IRS Custom Extractor (SPICE)| software to extract the spectra. 
With a flux limit of 5 mJy at 24 /im, we chose to use the optimal extraction with point-source 
calibration because it significantly improved the S/N ratios for our sources. When using the 
optimal method, each pixel was weighted by its position, based on the spatial profile of a 
bright calibration star. The outputs from SPICE produced one spectrum per order at each 
nod position, which were then combined. We also trimmed the ends of each order where 
the noise rose quickly. Finally, the flux-calibrated spectra of each order (including the 1st, 
2nd and 3rd orders) and module were merged without applying any scaling factor between 
SL and LL, and yielded a single spectrum per source. This spectrum was used to estimate 
aromatic feature fluxes, continuum flux densities at various wavelengths and line fluxes. 



2.3. Data Analysis 

2.3.1. The PAH Fluxes and Equivalent Widths 



To study the properties of PAH emission in our sample, we have used two methods to 
estimate the feature strength. The first method defines a local continuum or "plateau" under 
the emission features at 6.2 and 11.3 /xm by fitting a spline function to selected points, and 
measures the features above the continuum. The wavelength limits for the integration of the 
features are approximately 5.95-6.55 /im for the 6.2 /im PAH and 10.80-11.80 for the 11.3 /im 
PAH. We have not taken into account the possibility of water ice or HAC absorption in 
our measurement of the 6.2 /im PAH EW b ecause these featur es are known to be important 
mainly in strongly obscured local ULIRGs (ISpoon et al.l 120041 ); thus neglecting this compo- 
nent does not significantly change the 6.2 /im PAH EW. Although the 9.7/im silicate feature 
could affect the measurement on the 11.3/zm PAH, our sample has very few deeply obscured 
sources. The PAH EWs are derived by dividing the integrated flux over the average con- 
tinuum flux in each feature range. This PAH EW measured from the spline fitting method 
is defined as th e "apparent PAH EW" and is dire c tly comparable to th e stud i es in the lit- 
erature such as Peeters et al. fl2002h: ISpoon et all J2007h : krmus et al.l J2007h : besai et al.l 



software (jSmith et al. 



j2007h : IPope et al.l feoosh and bale et al.l J2009h . In the second method, we use the PAHFIT 



20071 ) to measure the PAHs in our sample (see Figured] for examples). 



1 For more details, check http://ssc.spitzer.caltech.edu/dataanalysistools/tools/irsclean/ 
2 For more details, check http://ssc.spitzer.caltech.edu/dataanalysistools/tools/spice/ 



-7- 



In PAHFIT, the PAH features are fit with Drude profiles, which have extend ed wings that 
account for a significant frac tion of the underlying pla teau ( jSmith et al.ll2007l ). As has been 
shown in lSmith et al.l ( 120071 ) and lGalliano et al.l ( 120081 ). although the PAHFIT method gives 
higher values of PAH integrated fluxes or EWs due to the lower continuum adopted than the 
"apparent PAH EW" method, the two methods yield consistent results on trends, such as 
the variations of band-to-band PAH luminosity ratios. Throughout this paper, when we refer 
to PAH EWs, we mean the apparent PAH EWs measured from the spline fitting method 
and they are used to classify object types. When we refer to PAH flux or luminosity, we 
mean the values derived from PAHFIT. 



2.3.2. The Fine- Structure Line Fluxes 

The mid-IR has a rich suite of fine-structure lines. [SIV]10.51/im, [NeII]12.81/zm, 
[NeIII]15.55/xm, [SIII] 18. 71/33. 48/zm and [SiII]34.82/xm are the most frequently detected 
fine-structure lines in the spectral range covered by the IRS. The high-excitation line of 
[OIV]25.89/im has often been detected in low metallicity galaxies, starburst galaxies or AGN, 
excited by the photoionization and/or shocks associated with intense star formation or nu- 
clear activity, while the [NeV] 14. 32/24. 32/zm lines are frequently detected in AGN-dominated 
sources and serve as unambiguous indicators of an AGN. 



We use the ISAP package in SMART (IHigdon et al.ll2004l ) to measure the strength of the 



fine-structure lines. A Gaussian profile is adopted to fit the lines above a local continuum. 
The continuum is derived by linear fitting except for the [Nell] 12. 81 /mi line, which is blended 
with the 12.7 //m PAH feature. The continuum underlying the [Nell] line is fit with a 2nd- 
order polynomial. The integrated fitted flux above the continuum is taken as the total flux 
of the line. Upper limits are derived by measuring the flux with a height of three times the 
local rms and a width equal to the instrument resolution. In this paper, we only use the flux 
ratio of [NeIII]/[NeII] to compare with the PAH strength, while the tabulated line fluxes will 
be presented and discussed in a future paper. 



3. The Infrared Luminosities of the 5MUSES Sample 



Several SED lib raries have been built to capture the variation in the shape of IR SEDs 
and t o estimate L m (bale fc Heloull2002l ; Icharv fc ElbazlhoOli braine fc Lill2007i : iRieke et al. 
20091 ). In the absence of multi- wavelength data, monochromatic luminosities have also 



been widely used to estimate L IR (jSajina et al.ll2007l ; iBavouzet et al.ll2008t IRieke et al.ll2009 
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Kartaltepe et al.ll2010l ). The 5MUSES sample has mid-IR spectra, in addition to the IRAC 



and MIPS photometry, which allows us to account properly for variations in the SED shape 
and obtain more accurate estimates of Ljr. 



3.1. Constructing an SED Template Library 



In order to cover a wide range of SED shapes to fit the 5MUSES sources, we have built 
an IR template library based on the recent observations obt ained from Spitzer. The library 
encompasses 83 ULIRGs observed by the IRS GTO sample (lArmus et al.l 120071); 75 normal 



star-fo rming galaxies from Spitzer Infrared N earby Galaxies Survey (SINGS. iKennicutt et al. 



(12003 )): and 136 PG and 2MASS quasars flShi et al.l 120071 ). The templates in the library 
consist of SEDs derived from IRS spectra and/or IRAC and MIPS photometry. For both 
the ULIRG and PG/2MASS sources, full 1-1000 /im SED have been obtained by Marshall 
et al. (2010, in preparation) and Shi e t al. (2010, i n prep aration) from IRS, MIPS and IRAS 
observations. For the SINGS galaxies, IDale et al.l (120071) have pr ovided SED fits to the MIPS 
24, 70 and 160 /im photometry using the IDale fc Heloul (120021 ) templates. However, these 
templates do not sample the full variation of the strength of PAH features in the 5-15 /im 
regime, due to the limited mid-IR spectra available when the templates were created. As 
a result, w h en we use SINGS galaxies as templates, we use their FIR SED from the fits of 
Dale et al.l ( 120071 ). while in the mid-IR, we use the observed IRAC photometry integrated 
from the whole galaxy. This extensive template library provides a good coverage on the 
variations of IR SEDs. 33% of our sources are best-fit with SINGS-type templates and 38% 
are best-fit with quasar-type templates. The remaining sources are best-fit by ULIRG-type 
templates. The type of the best-fit template also correlates well with the 6.2 /im PAH EWs. 
SB-dominated sources are normally best fit by SINGS-type templates and AGN-dominated 
sources are best-fit with quasar-type templates. For SB-AGN composite sources, the best-fit 
templates are divided among ULIRG, SINGS and quasar-type templates (48%, 37% and 
15% respectively). 



3.2. Estimating Ljr using Spitzer data 

Out of the 330 sources in 5MUSES, 280 galaxies have redshifts from optical or mid-IR 
spectroscopy. We are in the process of obtaining spec-z for the remaining 50 sources. We 
have estimated redshifts for 11 out of these 50 objects from silicate features or very weak 
PAH features, but do not include them in the discussion of this paper because of the large 
associated uncertainties. For the 280 objects, we use a combination of synthetic IRAC pho- 
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tometry obtained from the rest-frame IRS spectra, as well as the observed MIPS photometry 
to compare with the corresponding synthetic photometry from the SED templates and esti- 
mate total L IR . We select the best-fit template by minimizing \ 2 an d we use progressively 
more detailed and accurate Ltr estimation methods for 5MUSES source with more pho- 
tometry available. The final SED is composed of the IRS spectrum in the mid-IR and the 
best-fit template SED in the FIR. In the remainder of this section, we describe our method 
for estimating L IR and the associated uncertainties. 



3.2.1. Sources with MIPS FIR photometry 

For sources with FIR detection at MIPS 70 and 160 //m, we use five data points to 
fit their SEDs. The first two data points are rest-frame IRAC 5.8 and 8.0/iiiJ^ derived by 
convolving the rest-frame 5MUSES spectrum with the filter response curves of IRAC 5.8 and 
8.0 /im. The other three data points are the observed MIPS 24, 70 and 160 /im photometry 
for each 5MUSES source. The corresponding data points from the templates are derived 
in the following way: For ULIRG and PG/2MASS templates, the 5.8 and 8.0 /im fluxes are 
derived in the same manner as 5MUSES sources. The 24, 70 and 160 /im data points are 
derived by convolving the template SED at matching redshift with the MIPS 24, 70 and 
160 /im filter response curves. For SINGS templates, we use directly the observed IRAC 
5.8 and 8.0 /im photometry as the first two data points, which are essentially at rest-frame 
for all SINGS objects. Then we move the SINGS SEDs given by Dale et al. (2007) to the 
redshift of the 5MUSES source and derive the corresponding observed-frame MIPS 24, 70 and 
160 /im photometry. During the SED fitting, we weight the data points by their wavelength 
since the majority of the energy is emitted at FIR for IR selected sources, and look for the 
template that fits each 5MUSES source best by minimizing the \ 2 - A comparison of the 
ratio of flux densities at rest-frame 5.8, 8.0 /im and the observed MIPS 24, 70 and 160 /im 
photometry from the source and the best-fit template can be found in Figure [2] (solid line). 
The dispersion in the ratio of the observed photometry over the photometry from the best-fit 
template (F source /F temp i ate ) in each band is 0.07, 0.07, 0.03, 0.06 and 0.10 dex respectively. 

For sources with FIR detection only at MIPS 70 /im, we apply the same technique to fit 
the SED. We use the rest-frame IRAC 5.8 and 8.0 /im photometry and the observed MIPS 
24 and 70 /im data in our fitting. The upper limit at 160 /im is used to exclude templates 
for which the synthetic photometry exceeds the 2a upper limit of the source. A comparison 



3 5MUSES-312 has a redshift of 4.27 and for this source, we only use its MIPS 70 and 160 fim fluxes during 
the SED fitting. 
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of the ratio of flux densities at rest-frame 5.8, 8.0 /im and the observed MIPS 24 and 70 /im 
photometry from the source and the best-fit template can be found in Figure [2] (dashed line). 
The dispersion in the ratio of F source /F temp i ate in each band is 0.07, 0.07, 0.03 and 0.07 dex 
respectively. 

Once the best-fit template is identified, we derive Lir as explained in the last paragraph 
of the next section. 



3.2.2. Sources without MIPS FIR photometry 

19 objects do not have FIR detection even at 70 /im. For these sources, we select the IR 
SED based on the mid-IR spectra. Our method is to fit the IRS spectrum of the 5MUSES 
source with the mid-IR spectra of the templates in the corresponding wavelength regime and 
adopt the SED of the best-fit template. The templates of which the synthetic photometry 
exceed the 2a upper limits at MIPS 70 and 160 /im bands are excluded. As will be shown in 
Section 3.2.3, this IRS-only method might underestimate the Lir for cold sources by ~20%, 
while it shows no significant offset for warm sources^. All of the 19 objects in this category 
show SEDs with high (/W^m/iVoAtnOobs ratiog§; thus they are more likely to be warm sources. 
This suggests that our approach of using the IRS spectrum to find the best-fit SED is unlikely 
to result in significant biases on Ltr. 

Finally, for each source, we visually inspect the fitting results. We find that a range 
of templates could fit the SED well. We construct the final 5-1000 /im SED of a galaxy by 
combining its IRS spectrum in the mid-IR with the best-fit template SED at FIR. The total 
IR luminosity is derived by integrating under this SED curve. The uncertainty is derived 
from the standard deviation among the six best fits. We show examples of our SED fitting 
results in Figure [3] and the distribution of Lir is shown Figure H] a . We also show the 
distribution of L IR for each type of objects, e.g. starburst, AGN and composite (defined in 
detail in Section 4.1) in this figure. The derived Ltr for each source and its uncertainty is 
tabulated in Table [2J The distribution for the redshifts of 5MUSES objects are shown in 
Figure H]b. 



4 



Warm sources are d efined to have /24/jm/ '/Vo^m >0.2, derived from the definition of /25/xm//60Aim >0.2 



bv lSanders et all (|1988( ) 



5 



The 70 /Ltm flux densities are upper limits. 
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3.2.3. How well can we constrain IR SED from mid-IR? 



As has been shown above as well as in lKartaltepe et al.l ( 120101 ) . the availability of longer 



wavelength data greatly reduces the uncertainty in the estimate of Lir. We need to quantify 
how well one can constrain the SED of a galaxy if only the continuum shape up to ~30 /mi 
is available. In Figure El we show the comparison of the IRS predicted L{^ s and the L^ ot 
estimated from photometric data points (IRAC 5.8, 8.0/mi and MIPS 24, 70 and 160 /mi). 
For the IRS-only method, we use only the IRS spectrum and do not employ any longer 
wavelength information (70 and 160 /mi fluxes or upper limits) in our SED fitting, with the 
goal of testing solely the power of using mid-IR SED to predict FIR SED. We find that L IR 
estimated from mid-to-FIR photometry are on average 10% higher than Lir estimated from 
the IRS-only method, with a considerable scatter of 0.14 dex. It is worth noting that L{{^ s 
deviates from L^ ot by more than 0.2 dex for 20% of the sources while 5% of the sources 
deviate by more than 0.3 dex. We further divide the sources into two groups: cold sources 
and warm sources, based on the ratio of A^m/Ao^m- Cold sources (A^m/ Ao^m <0.2) 
show an average underestimate of 17% when using the IRS-only method and the la scatter 
is 0.16 dex, while warm sources (A^m/Ao^m >0.2) do not show systematic offset in the 
estimated Lir from the two methods, with a scatter of 0.18 dex. This comparison suggest 
that in IR-selected samples, the mid-IR spectrum could be used as an important indicator for 
Lir when no longer wavelength data are available. However, for cold sources, the IRS-only 
method might underestimate Lir by ~17%, due to the lack of information on the peak of the 
SED. For warm sources, although Lir estimated from the IRS-only method generally agrees 
with the value estimated from mid-to-far IR SED, the associated uncertainty is rather large. 
Thus for an individual galaxy, the L :R predicted by its mid-IR SED could be a factor of 1.5 
off from its intrinsic value for a significant subset of the population. 



3.3. Estimating the total IR Luminosity from a Single Band 

Using Spitzer data, we have obtained accurate estimates of the total infrared luminosi- 



often been used to estimate L TR ( 


3a una et al. 


2007; 


Papovich ei 


t al. 


2007; 


Pope et al. 


2008; 


Rieke et al. 


2009 




Bavouzet et al. 


2008 




Svmeonidis et al. 


2008 


). However, the fractional 



contribution of these photometric bands to the total infrared luminosity varies substantially 
depending on the dominant energy source. Because the IRS spectrum provides an unam- 
biguous way to identify the energy source for 5MUSES galaxies, our sample is ideal for 
investigating the difference in the fractional contributions of single band luminosities to Lir 
in different types of objects. 
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In Figure [61 we plot the ratio of several luminosity bands to Ltr. The PAH luminosities 
are plotted on the left panel and continuum luminosities are on the right. The dotted, dashed 
and dash-dotted lines respectively stand for the median ratios for SB, composite and AGN 
dominated sources. Clearly the fractional contribution of a certain band to Ltr is highly 
dependent on the type of the object, e.g. the monochromatic 24 fim continuum luminosity 
vh v accounts for ~13% of the total IR luminosity for SB galaxies, while it can contribute 
on average 30% of Ltr in AGN. The difference in the ratio of PAH luminosity to Ljr is less 
significant for different types of objects, because in order to be included on the left panel of 
this plot, the AGN dominated sources also need to have a solid detection of PAH feature 
that could be measured by PAHFIT, i.e. strong AGN sources are excluded. The mean ratios 
of Lgingieband/LiR are summarized in Table |3l Finally, we also provide our calibration of using 
single band luminosity to estimate Ljr in the Appendix. 



4. Aromatic Feature Diagnostics 
4.1. The Average Spectra 



We derive the stacked SEDs for the SB, composite and AGN dominated sources in the 
5MUSES sample, combining the low resolution IRS spectra in the mid-IR with the MIPS 
photometry at FIR. Alth ough we do not have optical spectroscopy to classify the object types 
with the BPT diagram ( [Baldwin et al.lll98ll ; iKewley et al.ll200ll ). the equivalent widths of 
PAH features can be used as indicators of star formation activity. The 6.2 and 11.3/im 
PAH bands are relatively isolated with little contamination from nearby features, which is 
important for unambiguously defining the local continuum. However, the 11.3 /xm band is 
located on the shoulder of 9.7/im silicate feature. Thus its integrated flux and underlying 
continuum are likely to be affected by dust extinction effects. As a result, in our discussion, 
we use the 6.2 /im PAH EWs to classify objects. To be consistent with the studies in the 
literature, we have adopted the following criteria for our spectral classification: sources with 
EW>0.5 /im are SB-dominated; sources with 0.2< EW<0.5 /im are A GN-SB composite and 
sources with EWs < 0.2 /jm are AGN-dominatecfl. ( Armus et aDl2007l ). The PAH EWs for 
the sample are tabulated in Table 121 Out of the 280 sources for which redshifts have been 
obtained from optical or infrared spectroscopy, there are 123 SB galaxies (44%), 62 composite 
sources (22%) and 95 AGN dominated sources (34%). 



6 Sources with a significant old stellar population could also have a reduced 6.2 /im PAH EW. As will 
be shown in Shi et al. 2010 (in preparation), the stellar emission contributes less than ~20% to the 6/im 
continuum for our IR selected sample of 5MUSES. 
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The 5-30 /mi composite spectra are derived by first normalizing individual spectra at 
rest-frame 5.8 /mi, and then taking the median in each wavelength bin. In Figure [7J a, we 
show the typical SED for SB galaxy in blue and AGN in red, while yellow line represents the 
median SED for SB- AGN composite sources in 5MUSES. The average SEDs have been offset 
vertically. The shaded regions represent the 16th and 84th percentile of the flux densities at 
each wavelength. 

The MIPS 70 and 160 /mi photometry is crucial for constraining the SED shape of 
a galaxy and we have also included these data in the final typical SED. Because of the 
difference in redshift range for sources of different spectral types, we have divided the MIPS 
70 and 160 /mi data into several rest-wavelength bins before we take the median. For SB and 
composite sources, we take 2 bins: 40-70 /mi and 70-160 /mi; For AGN, we choose to have 3 
bins due to their larger redshift range: 30-50, 50-100, and 100-160 /mi. We take the median 
flux in each bin and assign the 16th and 84th percentile of the data points in the same bin 
as the uncertainties. The final median SEDs are presented in Figure UJb. We can clearly see 
that besides having much less PAH emission in the mid-IR, the continuum in the AGN also 
rises much more slowly than in the SB. The SED of the composite source is between the SB 
and AGN and its shape is dependent on how we define a composite source. As can be seen 
in Figure [7] and EJ our definition of composite sources with 0.2 /mi<6.2 /mi PAH EW<0.5 /mi 
is likely biased towards star formation dominated sources (see Section 4.2). 



4.2. The Distribution of PAH EWs 

With the superb sensitivity and spectral coverage of the IRS, we are able to quantify the 
strength of the PAH emission over nearly two orders of magnitude in its EW. The distribution 
of the 6.2 /mi PAH EWs for the 280 known-redshift galaxies in 5MUSES is shown in Figure 
[HI The solid line represents the distribution for sources with detection of the 6.2 /mi feature, 
while the dotted line also includes upper limits. We clearly observe a bimodal distribution in 
Figure El with two local peaks at ~0.1 and ~0.6/xm. This is somewhat surprising, because 
5MUSES provides a representative sample completely selected based on IR flux densities, 
and one would have expected a more continuous distribution. Although we still lack redshift 
information for 50 sources in our sample, the featureless power-law shape of their IRS spectra 
(except for a few cases where silicate absorption or very weak PAH feature is present) indicate 
that these are likely to be AGN-dominated. Thus if they were included in Figure [HI they 
would most likely be located in the range between 0—0.2 /mi, and the bimodal distribution 
would not be affected. A similar bi-modality is also observed in the distribution of the 
11.3 /mi PAH EWs (not shown here). The observed bimodal distribution of the PAH EWs 



-14- 



may be a result of the selection effect for this flux limited sample: objects at higher redshifts 
are more likely to be AGN and thus pile up at the low EW end. However, if we divide our 
sample into sources with z>0.5 and z<0.5, the bimodal distribution is again observed in the 
z<0.5 population, although all the z>0.5 objects are located at the low EW end. Detailed 
population modelling is being performed and this issue will be addressed in a later paper. 



4.3. PAH Properties versus mid-IR and FIR slopes 



Another important physical parameter that is often used to quantify the dominant 
energy s ource of a galaxy is the ratio of w arm to cold dust. It has been shown in previous 
studies JDesai et al.ll2007t IWu et aI.ll20oih that the 6.2 and 11.3 fim PAH EWs of galaxies 
are usually suppressed in warmer systems dominated by AGN, as indicated by the low flux 
ratios of IRAS fm/f 25- For the 5MUSES sample, we have examined the correlation of the 
6.2 jum PAH EWs with various continuum slopes, e.g. fis/fs.s, ho/h.8, ho/hs and ijo/^i- 
The rest-frame continuum fluxes are estimated from the final SED obtained from the fits in 
Section 3. We find that the continuum ratios of f3o/fis and f7o/f24 have the strongest global 
correlation with the 6.2 /zm PAH EWs and the correlation coefficients are both ~0.7. 

In Figure IH] a and b, we plot the 6.2 /im PAH EW against /30//15 and fyo/izi- The 
5MUSES populations separate into two groups, one with steep spectra and high aromatic 
content, and the other with slow rising spectra and low aromatic content. The gap between 
SB and AGN dominated sources is likely due to the selection effect of this sample. We should 
note that within each group, there is little if any correlation between the slope and the PAH 
EW, but it is the contrast between the two grou ps that gives t h e ove rall impression of a 
correlation. This is consistent with the studies of IVeilleux et all (120091 ). who have showed 
the power of using the 7.7 //m PAH EWs and /30//15 ratios as indicators of AGN activity, 
despite the large scatter associated with each parameter. To understand the variation in the 
PAH EWs and continuum slopes, we further divide our sample into smaller bins and estimate 
the average values in each bin. The sources are divided according to their f7o/f24 ratios or 
f3o/fi5 ratios and we assign an equal number of objects to each bin. We find that sources 
in the first three bins with log(f3o/fis)>0.65 or Iog(f7o/f24)> 0.73 all have median 6.2/im 
PAH EWs of ~0.60/zm and dispersion of ~0.2dex, which again confirms our observation 
that within the group of starburst galaxies, there is little correlation between the slope 



7 If the spectral index is defined as a=log(fl/f2)/log(Vl/^2), then the continuum slope ratios of 
Iog(i3o/fi5)>0.65 can be translated to a 30 -i5 <-2.17 and log(fYo/f24)> 0.73 can be converted to 070-24 <- 
1.57. 
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and aromatic content. Sources with Iog(f3o/fi5)<0.38 or Iog(f7 /f24)<0.38 are clearly AGN- 
dominated with very low PAH EW. The median values and the associated uncertainties of 
the 6.2 fim PAH EWs and continuum ratios are summarized in Table HI 

We also investigate the variation in the ratio of Lpah/Ltr when the galaxy color indi- 
cated by the continuum slope changes. We use the sum of the PAH luminosity from the 
6.2 fim, 7.7 (Urn complex and 11.3 fim complex to represent Lpah ; measured from the compos- 
ite spectra derived for each bin. It is clear that the PAH fraction stays nearly constant for 
starburst dominated systems, while its contribution drops significantly when AGN becomes 



more dominant (see Table H]) . 



t has been shown that the PAH luminosity can contribute 



10% in star-forming galaxies (ISmith et al.l 120071 ). For our sample, we find that Lpah con- 



tributes ~5% to L IR . This ratio is lower than the SINGS results. We have only taken the 6.2, 
7.7 and 11.3 /iin bands into account^], while the SINGS studies include all the PAH emitting 
bands in the mid-I R. Since the 6.2, 7 .7 and 11.3 /im bands accounts for ~68% of the power 



in PAH emission ( ISmith et al.l 120071 ). our Lpah/Lir ratio can be converted to ~7.5% for 
the total PAH contribution to Lpp. This is still slightly lower than the SINGS results, but 
consistent within uncertainties. 

Finally, for each group of continuum slope sorted spectra, we derive typical 5-30 fim 
SEDs by taking the median flux densities in every wavelength bin after normalizing at rest- 
frame 5.8 fim. This will be useful for SED studies when only galaxy colors estimated from 
broad band photometry are available. These composite SEDs are shown in Figure [TDJ Then 
we explore whether the derivation of total IR luminosity from broadband photometry varies 
with galaxy color. We assume Lip is correlated with L 2 4 /im and L 70/jm in the following manner 
and derive the a and b coefficients in each f7o/f24 continuum slope bin (all in rest-frame): 

logL IR = a \ogL 2ifim + b \ogL 70fim (1) 

The values of a and b coefficients are summarized in Table HI To illustrate the variations 
in each slope bin, we plot the ratio of Lip/L 2 4 Mm versus Lyo^m/L^m in Figure [11] a. The 
sources are colored according to their fyo/f24 ratios. We clearly observe that when normalized 
by the monochromatic 24 /im luminosity, Ljr is strongly correlated with L 70 ^ m and the slopes 
in each continuum ratio bin become steeper when L 70Mm /L 2 4 Atm increase, except in the last 
slope bin (see also the b coefficients). We fit a 2nd-order polynomial to the data and find 
the correlation to be: 

logy^- = (0.476 ± 0.005) + (0.509 ± 0.010) log^^ + (0.370 ± 0.022) (log^^) 2 (2) 

-^24/xm -^24/^m 



8 The S/N ratios of the 5MUSES spectra are much lower than SINGS, thus we only include the strongest 
PAH bands. 
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The above equation is derived based on the 5MUSES data. The majority (90%) of the 
24 yum luminosities of these 280 galaxies are between 10 9 L Q and 10 12 L Q . The ratio of 
f7n / t m/f?4/Am ranges fr o m 0.4 5 to 34. Our result is consistent with a similar correlation derived 



by iPapovich &: Belli (120021 ). while it diverges for sources with low ^o^/h^ra ratios, since 
their modelling work has focused on star-forming galaxies only. 

We repeat the same exercise for our sample binned with the fso/fis ratios. In Figure [TT1 
b, we find that L IR is correlated with L 30Atm when both quantities are normalized by L 15/mi , 
although with very large scatter. The dotted line is a linear fit to the data. For a given 
L30/im/Li5/im ratio, LiR,/Li5 Mm can span as much as a factor of five. The median values in 
each group binned by the fso/fi5 ratios are also summarized in Table HI 



4.4. The Variation in PAH Band-to-Band Strength Ratios 

The luminosity ratio of different PA H bands is thought to be a function of the grain 
size and ionization state ( Tielens 2008 ). luminosity ratios of LpAH7.7^m/LpAHii.3^nS with 
the 6.2 /im PAH EWs for the 5MUSES sample. Only sources with S/N>3 from PAHFIT 
measurements for the 7.7 and 11.3 /im bands are included in this plot. We find that the AGN- 
dominated sources on average have lower L P AH7.7^m/LpAHii.3/xm ratios than the composite 
or SB-dominated sources. The mean log(L P AH7.7^m/LpAHii.3/xm) ratios for AGN, composite 
and SB galaxies in 5MUSES are 0.32±0.18, 0.53±0.15 and 0.53±0.08, respectively. This is 
consistent wit h the studies on the nuclear spectra of low luminosity star-forming galaxies 



from SINGS (ISmith et al.l 120071). which a ls o sho w decreased L P AH7.7^m/LpAHii.3^m ratios in 



spectra with AGN signals. ISmith et al.l (120071 ) suggest that this change in the ratio of 
LpAH7.7^m/LpAHii.3^m is likely due to the destruction of the smallest PAHs by hard photons 
from the AGN. On the other hand, AGN are less extinguished than SB or composite sources, 
thus if PAHFIT underestimates the extinction correction, it will preferentially underestimate 
the 11.3 /xm fluxes more than the 7.7 jum feature in SB/Composite sources, thus resulting in 
the elevated ratios of L P AH7.7/xm/LpAHii.3^m m SB/composite systems. 



In Figure [13] a, we show the histogram of the L P AH7.7/xm/L P AHii.3/xm ratios for the SB- 
dominated sources in 5MUSES. We have overplotted the values from the SINGS sample. 
To make a fair comparison, we remeasure the PAH luminosity and EWs for the SINGS 
nuclear spectra using the same method as 5MUSES and classify the sources with 6.2 /iin 
PAH EWs larger than 0.5 /im as SB-dominated. We have also included the distribution 



9 We choose to use the Lpah7 7L tm/LpAHii.3um ratio in this study for easier comparison with literature 
results, such as lSmith et all (|2007t l: lO'Dowd et al.l |2009l ). 
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of the LpAH7.7//r T l /LpAHii.,ym ratios f rom the UV/SDSS selected star-forming galaxies sam- 
ple of SSGSS (jO'Dowd et all 120091 ) . For this last sample, the star- forming galaxies are 
classified from opt ical spectroscopy using the BPT diagram method ([Baldwin et all Il98ll : 
Kewlev et aJboOlh . We find that the distribution for SB galaxies in 5MUSES and SSGSS 
is similar, while both samples appear to have lower LpAH7.7/mi/LpAHii.3/mi ratios than the 
nuclear spectra of SINGS SB galaxies. The mean log(LpAH7.7^m/LpAHii.3/mi) ratio for SINGS 
starbursts is 0.63±0.06 while it is 0.53±0.08 for 5MUSES starbursts. This might be a reso- 
lution effect: If the physical conditions at the nuclear region of a galaxy indeed modifies the 
distribution of the LpAH7.7/mi/LpAHii.3/mi ratios, it might be visible only in the spectra taken 
through apertures with small projected sizes. The median redshift for the SB-dominated 
sources in the 5MUSES sample is 0.12, while the median redshift for the SSGSS sample is 
0.08. At the redshift of 0.08, 1" corresponds to 1.53 kpc. The IRS spectra (for SL, the slit 
width is~3.6") of 5MUSES and SSGSS source s are integra t ed fro m the whole galaxy, thus 
diluting the signature of the nuclear regions. ISmith et all ( 120071 ) have shown the changes 
in the LpAH7.7^m/LpAHii.3^m ratios in spectra extracted from bigger to smaller apertures 
in two star-forming galaxies: The LpAH7.7/ tm /LpAHii.3^m ratios measured from star-forming 
galaxy spectra extracted with smaller apertures are hi gher than those measured from larger 



apertures, consistent with our results. More recently, iPereira-Santaella et al. I (120 101 ) have 
suggested that the 11.3 /zm PAH feature is more extended than the 6.2 or 7.7/zm PAH from 
a spatially resolved mapping study of local luminous infrared galaxies. They have observed 
lower LpAH6.2/im/LpAHii.3/xm ratios in the nucleus, consistent with our results. We also show 
the distribution of LpAH6.2/mi/LpAH7.7/xm ratios in Figure [T3] b. No significant difference has 
been observed between the 5MUSES, SINGS and SSGSS samples. 



Finally in Figure [HJ we present the variation in PAH band-to-band ratios for the three 
strongest bands at 6.2, 7.7 and 11.3/zm of the 5MUSES sample. Only sources with S/N>3 
from PAHFIT measurements for all three PAH bands are included in this figure. The 
two dark lines represent the traces for fully neutral or fully ioniz ed PAH molecules with 
different numbers of carbon atoms predicted from modelling work (jDraine fc Lill200lf ). The 
LpAH7.7|tm/LpAHii.3|tm ratios span a range of a factor of 5 while the LpAH6.2/Mi/LpAH7.7,um ratios 
only vary by a factor of 2. The uncertainty in the L P AH7.7^m/LpAHii.3/xm ratios is 0.09 dex and 
it is 0.05 dex for the Lp A H6.2^m/LpAH7.7 M m ratios. This narrow ran ge of Lpahr.^/Lp ah? 7,, m 
ratios is consistent with the values for the SINGS nuclear sample (ISmith et al.ll2007l ). while 
we have not observed any sources wit h extremely low LpA Hfi.2/i.m/LpAH7.7/mi ratios (<0.2) as 
has been found in the SSGSS sample (jO'Dowd et al.ll2009l ). 
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4.5. PAH band ratio versus [NeIII]/[NeII] 



Because of the large difference in ionization potentials of the Ne ++ (41eV) and Ne + 
(21.6eV) ions, the ratio of [NeIII]/[NeII] is often used as a tracer of the hardness of the 
radiation field. The [Nelll] 15.55 /zm and [Nell] 12.81 fim lines are among the strongest lines 
emitted in the mid-IR and because differential extinction effects between their wavelengths 
are small, they are particularly valuable. We use the IRS low-res spectra to identify and 
measure these lined_|. The line fluxes measured from low-resolution spectra have on average 
an uncertainty of ~20%. 

In Figure [T5J we show the flux ratios of L P AH7.7^m/LpAHii.3^m versus [NeIII]/[NeII]. 
The solid symbols denote detections while the open triangles represent upper/lower lim- 
its. We overplot the median LpAH7.7^m/LpAHii.3^m for SB-dominated sources in 5MUSES as 
the dotted line. We find that the SB, composite and AGN-dominated sources (including 
sources with upper/lower limits) are almost evenly distributed on the two sides of the dot- 
ted line. However, the AGN with solid detections on both axes do appear to have lower 
LpAH7.7/im/LpAHii.3^m ra tios in gener a l. As has been discussed in Section 4.4, this is consis- 
tent with the studies of ISmith et all ( 120071 ) using the SINGS nuclear spectra. We note that 
the 5MUSES sample do not have sources with extreme LpAH7.7^m/LpAmi.3/xm ratios compa- 
rable to the lowest ones reached by SINGS. This is probably because the SINGS spectra 
probe smaller, more central and thus more AGN-dominated regions. It should also be noted 
that the AGN lum inosities in 5MUSES are substantially higher than SINGS. Our results 
are consistent with lO'Dowd et all (120091 ). who have studied a UV-SDSS selected sample at 
z~0.1 and do not observe extreme LpAH7.7/mi/LpAHii.3/mi ratios either. We also notice that the 
range of [NeIII]/[NeII] ratios are similar for all three groups of objects t hat we have classified 



based on their 6.2 /mi PAH EWs. This is consistent with the study of iBernard-Salas et al. 



(120091 ). who found no correlation between the PAH EWs and the [NeIII]/[NeII] ratios in 
a sample of starburst galaxies. However, in more extreme radiation field conditions, such 
as low-metallicity environment, PAH EWs have been ob served to anti-c orrelate with the 
radiation field hardness indicated by [NeIII]/[NeII] ratios (jWu et all 120061 ) . 



10 For the 5MUSES sample, only 21 out of 330 sources have IRS high-resolution spectra, which limits our 
ability to probe the full dynamic range covered by the whole sample. Thus we use the low-resolution spectra 
to measure the [Nell] and [Nelll] fluxes to compare with the PAH band-to-band ratios. 
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5. Conclusions 

We have studied a flux limited (i^nm >5mJy) representative sample of 330 galaxies 
surveyed with the Infrared Spectrograph on board the Spitzer Space Telescope. Secure 
redshifts of 280 objects have been obtained from optical or infrared spectroscopy. The 
redshifts of the 5MUSES sample ranges from 0.08 to 4.27, with a median value of 0.144. 
This places the 5MUSES sample at intermediate redshift, which bridges the gap between 
the nearby bright sources known from previous studies and the z~2 objects pursued in 
most of the IRS follow up observations of deep 24 /zm surveys. The simple selection criteria 
ensures that our sample provides a complete census of galaxies with crucial information on 
understanding the galaxy evolution processes. 

Using mid-IR spectroscopy and mid-to-far IR photometry, we have obtained accurate 
estimates on the total infrared luminosities of 5MUSES galaxies. This is achieved by min- 
imizing the x 2 to find the best fit template from our newly constructed empirical SED 
library built upon recent Spitzer observations. The availability of longer wavelength data 
also greatly reduces the uncertainties in Lip. When only one IRS spectrum is available, 
one can still predict the shape of the FIR SED from the mid-IR and estimate Lir, albeit 
with substantially larger uncertainties (0.2 dex). The IRS-only method does not introduce a 
systematic bias when estimating L IR for warm sources, but could underestimate the L IR by 
~17% for cold sources, due to the lack of information sampling the peak of the SED. The 
fractional contribution of single band luminosity to Lir varies depending on the dominant 
energy source and the average values have been calculated for the SB, composite and AGN 
dominated sources, as well as the whole sample. 

We analyze the properties of the PAH emission in our sample using the IRS spectra. 
The PAH EWs show a bimodal distribution, which might be related to the selection effect of 
the sample. The starburst and AGN dominated sources form two clumps when comparing 
the continuum slopes and PAH EWs, while there is little discernible correlation within 
each group. Average spectra binned with the 6.2 /im PAH EWs, the continuum slopes of 
log(f 30 /fi5) and Iog(f7o/f24) have been derived to show the typical SED shapes. The variation 
in PAH EW and Lpah/Lfr ratios when galaxy color changes have also been inspected. The 
galaxy color provides essential constraint on estimating the total infrared luminosity from 
broadband photometry. 

We have also inspected the band-to-band PAH intensity ratios with regard to different 
spectral types. The LpAH7.7 M m/LpAHii.3 M m ratios in AGN dominated sources in 5MUSES are 
on average lower than the SB or composite sources. The SB, composite and AGN dominated 
sources have mean log(L PA H7.7^m/LpAHii.3^m) ratios of 0.53±0.08, 0.54±0.15 and 0.32±0.18, 
respectively. The mean log(L P AH7.7//m/LpAHii.3/im) ratio for the SB dominated sources in 
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5MUSES is lower than the mean ratio derived from the nuclear spectra of SB galaxies in 
SINGS (0.63±0.06), which might indicate a difference in the physical conditions near the 
nucleus versus over the entire galaxy. At the median redshift of our sample, the IRS SL slit 
width corresponds to a few kpc, thus even if the ionization state or grain size distribution 
is different at the nuclear level, the signal might get diluted when we study the integrated 
spectrum and would result in the different log(LpAH7.7 ( um/LpAHii.3 / um) ratio distribution. 

Finally, we provide our calibration of using PAH luminosity or mid-IR continuum lumi- 
nosity to estimate Ljr in the Appendix. We have shown that single band luminosities trace 
the Lir differently in SB or AGN dominated sources and we provide calibrations for each 
object type. This technique will be useful for luminosity estimates when no mult i- wavelength 
data are available. 

We thank the anonymous referee whose comments have helped to improve this manuscript. 
This work was based on observations made with the Spitzer Space Telescope, which is op- 
erated by JPL/Caltech under a contract with NASA. The observations are associated with 
the Spitzer Legacy Program 40539. The authors acknowledge support by NASA through 
awards issued by JPL/Caltech. This research has made use of the NASA/IPAC Extragalac- 
tic Database (NED) which is operated by the Jet Propulsion Laboratory, California Institute 
of Technology, under contract with the National Aeronautics and Space Administration. 

A. Estimating the Total Infrared Luminosity from PAH or Monochromatic 

Continuum Luminosities 

In Section 3, we have discussed in detail our method to estimate the total infrared 
luminosities for the 5MUSES sources. The empirical library of SED templates built from 
Spitzer observations, as well as the availability of photometric and spectroscopic data from 
mid-IR to FIR for 5MUSES, allow us to have precise estimates on their Lir. We have shown 
in Figure E] the importance of having FIR data in determining the total energy output in 
the infrared. However, for high redshift galaxies, FIR observations are not always available. 
Herschel Space Observatory will provide FIR measurements from 70 to 500 /im to reveal the 
properties of cold dust in many systems. For now, we provide our calibration of estimating 
Lir from several bands in the mid-IR and discuss its applications. The following correlations 
are derived by performing a linear fit to the 5MUSES data with equal weight on each object 
because the dispersion of the data point in the x-y plane is larger than the measurement 
errors. 

As has been shown in many studies, the infrared SED of a starburst galaxy is drastically 
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different from that of an AGN flBrandl et all 12009 : lHao et al.l 120071 : lArmus et all 120071 ). 



Because of these substantial variations in the SED shapes, it is crucial to calibrate the 
luminosity estimates for each spectral type. Here we provide our luminosity calibrations 
based on the three spectral types : starburst, composite and AGN. The following PAH 
luminosities are derived from the PAHFIT method. 

1. 6.2 jum PAH: With a wavelength cut at 28 /im for the James Webb Space Telescope 
(JWST), the 6.2 PAH feature might be the only PAH band that could be observed to 
quantify star formation activities in z~3 sources when JWST is launched. 

For SB sources: 

logL IR = (2.40 ± 0.22) + (0.96 ± 0.03)logL PAH6 . 2/ ,m (Al) 
For composite sources: 

logL IR = (1.76 ± 0.32) + (1.04 ± 0.04)logL PAH6 . 2/ , m (A2) 

For AGN sources: 

logL IR = (-0.58 ± 0.58) + (1.30 ± 0.06)logL PAH6 .2 M m (A3) 

2. 7.7 jum PAH: The 7.7 /im PAH complex is the strongest band among the various PAH 
features. It is often used to estimate the total infrared luminosities for the z~l-2 sources 
pursued in IRS observations of 24 /im selected sources. 

For SB sources: 

logL IR = (2.30 ± 0.26) + (0.91 ± 0.003)logL PAH7 . 7Atm (A4) 
For composite sources: 

logL IR = (1.80 ± 0.44) + (0.98 ± 0.05)logL PAH7 ^m (A5) 

For AGN sources: 

logL IR = (3.45 ± 0.90) + (0.83 ± 0.09)logL PAH7 . 7 ^ m (A6) 

3. 11.3 /im PAH: The 11.3 /im band is another strong PAH band in the mid-IR that 
is relatively isolated from other PAH bands. However, the integrated fluxes from this band 
might be affected by the 9.7 /xm silicate feature. 

For SB sources: 

logL IR = (2.18 ± 0.22) + (0.98 ± 0.03)logL PAH n. 3 ^ m (A7) 
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For composite sources: 

logL IR = (1.49 ± 0.43) + (1.07 ± 0.05)logL PAH ii.3 M m 

For AGN sources: 



logL m = (2.22 ± 0.36) + (1.00 ± 0.04)logL PA Hii 



(A? 



(A9) 



4. 6.2+7.7+11.3 /im PAH: In normal star-for ming galaxies, the PAH emission accounts 
for ~10%-15% of the tota l infrared luminosi ties ( jSmith et al.l 120071 ). while this fraction is 
smaller for local ULIRGs (lArmus et al.l 120071 ). Here we use the sum of the three strongest 
PAH bands, the 6.2, 7.7 and 11.3 /xm PAH luminosities to represent the total PAH lumi- 
nosities. However, when using the correlation provided here, one needs to keep in mind that 
the properties of PAHs studied in the local universe might be different at hi gh z, as has 



already been revealed in the study of several z~2 luminous infrared galaxies (ISajina et al. 



20071 : iPope et al.l 120081 ) . Understanding the PAH contribution in our intermediate redshift 
sample would also be instrumental for tackling the problem of whether and how PAH emis- 
sion evolves with redshift in future studies. 
For SB sources: 



logL IR = (1.82 ± 0.25) + (0.95 ± 0.03)logL PAH 6.2 + 7.7+ii.3 Mm 
For composite sources: 

logL IR = (0.73 ± 0.36) + (1.06 ± 0.04)logLp A H6.2+7.7+ii.3 / ,m 
For AGN sources: 



logL IR = (-2.13 ± 1.01) + (1.35 ± 0.10)logLpAH6.2 +7 .7+i 



1.3/xm 



(A10) 



(All) 



(A12) 



5. 5.8 fim monochromatic continuum luminosity: The 5.8 fim continuum luminosity 
provides a crude estimate of Li R . In AGN dominated sources, the 5.8 fim continuum will 
be elevated due to the presence of very hot dust component. This is also a band that is 
available for most of the high redshift samples observed by Spitzer, and for JWST when it 
is launched. 
For SB sources: 

logL IR = (1.94 ± 0.25) + (0.95 ± 0.03)logL 5 . 8/tm (A13) 

For composite sources: 



logL IR = (2.68 ± 0.45) + (0.87 ± 0.05)logL 



5.8/jm 



(A14) 
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For AGN sources: 

logL m = (2.34 ± 0.29) + (0.85 ± 0.03)logL 5 . 8Mm (A15) 

6. IRAC 8 jum: The rest-frame IRAC 8.0 /iin band has included both dust continuum 
emission and PAH emission from the 7.7, 8.3 and 8.6 /xm PAH band (if present). It provides 
a useful channel for estimating Lir from PAH features when no spectroscopy is available. 
For SB sources: 

logL IR = (1.60 ± 0.21) + (0.93 ± 0.02)logL IRAC8 ^ m (A16) 
For composite sources: 

logL IR = (1.45 ± 0.26) + (0.95 ± 0.03)logL I R AC 8^m (A17) 

For AGN sources: 

logL IR = (1.70 ± 0.26) + (0.90 ± 0.02)logL I R AC 8^m (A18) 

7. 14 /im monochromatic continuum luminosity: The 14 /iin is an important band in 
the mid-IR that is still sensitive to the AGN emission. 

For SB sources: 

logL IR = (1.44 ± 0.16) + (0.97 ± 0.02)logL 14/im (A19) 

For composite sources: 

logL IR = (1.76 ± 0.37) + (0.93 ± 0.04)logL 14/im (A20) 

For AGN sources: 

logLiR = (1.62 ± 0.29) + (0.90 ± 0.03)logL 14 ^ m (A21) 

8. 24 /im monochromatic continuum luminosity: Here we refer to the 24 /im continuum 
luminosity averaged in one micron range, instead of the rest-frame MIPS 24 /im band. This is 
because if we use the MIPS 24 /im band, sources at z>~0.3 will be eliminated from this study 
due to the limited wavelength coverage of its rest-frame mid-IR spectra. The sources we use 
in the calibration mostly have 10 10 L o <L IR <10 12 L Q and no quasars have been included in 
this calibration because of the wavelength cut. Since our sample is selected at 24 /im, it tends 
to favor warmer sources, which also needs to be kept in mind when using these relations. 
For SB sources: 

logL IR = (1.56 ± 0.23) + (0.93 ± 0.02)logL 24/im (A22) 
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For composite sources: 

logL IR = (1.64 ± 0.44) + (0.91 ± 0.04)logL 24/im (A23) 

For AGN sources: 

logL IR = (1.67 ± 0.47) + (0.89 ± 0.05)logL 24 ^ m (A24) 
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Table 1. On-Source Integration Time of the Sample 



f24/im ( m Jy) SL2 (second) SL1 (second) LL2 (second) LLl (second) 



5~7 480 480 480 480 

7-10 480 240 240 240 

10~15 480 480 180 180 

15~25 240 120 120 120 

>25 120 120 60 60 
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Table 2. General Properties of the Sample 



ID Name RA (J2000) Dec (J2000) Rcdshift a f 2 4 M m (mjy) 6.2^m EW log(L m /L ) 
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05 


5MUSES 


319 


5MUSES 


.J172159 


43+595034 


3 


17h21m59 


3s 


+59d50m34 


2s 





028(2) 


9 


7 


0.387+0 


031 


9 


78±0 


03 


5MUSES- 


320 


5MUSES 


.J172219 


58+594506 


9 


17h22ml9 


6s 


+59d45m07 


Os 





272(2) 


7 


8 


<0.133 




11 


24±0 


02 


5MUSES 


321 


5MUSES 


.J172228 


04+601526 





17h22m28 


2s 


+60dl5m26 


2s 





742(2) 


7 


2 


<0.111 




12 


40±0 


03 


5MUSES- 


322 


5MUSES 


.J172238 


73+585107 





17h22m38 


8s 


+58d51m07 


Os 


1 


624(1) 


6 


7 


<0.062 




13 


12±0 


04 


5MUSES- 


323 


5MUSES 


.J172313 


06+590533 


1 


17h23ml3 


Is 


+59d05m33 


Is 





108(2) 


6 


2 


0.750+0 


037 


10 


85±0 


03 


5MUSES- 


324 


5MUSES 


J172355 


58+601301 


7 


17h23m55 


5s 


+60dl3m01 


Is 





175(2) 


5 


4 


0.905+0 


034 


11 


13±0 


02 


5MUSES 


325 


5MUSES 


.J172355 


97+594047 


6 


17h23m56 


Os 


+59d40m47 


4s 





030(2) 


5 


2 


0.518+0 


098 


9 


35±0 


04 
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Table 2 — Continued 



ID 


Name 


RA (J2000) 


Dec (J2000) 


Redshift a 


hi^m (mjy) 


6.2^m EW 


log(L ra /L ) 


5MUSES-326 
5MUSES-328 
5MUSES-329 
5MUSES-330 


5MUSES.J172402. 1 1+600601 .4 
5MUSES.J172546.80+593655.3 
5MUSES.J172551. 35+601138.9 
5MUSES.J172619.80+601600. 1 


17h24m02.1s 
17h25m46.8s 
17h25m51.3s 
17h26ml9.8s 


+60d06m01.2s 
+59d36m55.3s 
+60dllm38.9s 
+60dl6m00.0s 


0.156(2) 
0.035(2) 
0.029(1) 
0.924(1) 


8.0 
26.0 
27.3 

6.5 


0.461+0.024 
0.554+0.041 
0.454+0.005 
<0.039 


11.13+0.03 
10.49+0.04 
10.25+0.03 
12.35+0.08 



a Thc rcdshifts obtained from NASA/IPAC Extragalactic Database arc indicated with "1" while the redshifts derived from IRS spectra 
arc indicated with "2" . 



Table 3. Median Luminosity Ratios of the Sample 





SB 


composite 


AGN 


Whole Sample 


i°g(LpAH6.2/tm/LlR) 


-2.03±0.13 a 


-2.11±0.13 


-2.22±0.14 


-2.06±0.14 


l0g(LpAH7.7/im/LlR,) 


-1.51±0.15 


-1.54±0.18 


-1.76±0.28 


-1.53±0.20 


log(LpAH11.3,ttm/LlR) 


-2.04±0.13 


-2.14±0.20 


-2.20±0.18 


-2.08±0.17 


log(LpAH6.2+7.7+11.3,um/LlpJ 


-1.29±0.14 


-1.33±0.12 


-1.39±0.19 


-1.31±0.14 


log(L 5 . 8Mm /L IR ) 


-1.44±0.15 


-1.36±0.29 


-0.77±0.32 


-1.33±0.42 


log(Li R AC8Atm/LlR) 


-0.90±0.13 


-0.96±0.15 


-0.66±0.23 


-0.86±0.22 


log(L 14/im /L IR ) 


-1.19±0.10 


-1.11±0.21 


-0.57±0.19 


-1.12±0.31 


log(L 24M m/LlR) 


-0.87±0.14 


-0.82±0.23 


-0.49±0.15 


-0.81±0.22 



a The dispersion is the la deviation for each group of objects. 
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Table 4. Median PAH strengths and Continuum Ratios of the Sample 



log(f 3 o/fi 5 ) a 


log(f 70 /f 2 4) 


6.2 ^m PAH EW b 


log(LpAH/Ll R ) C 


a 


b 


>0.898 
0.793-0.898 
0.650-0.793 
0.376-0.650 

<0.376 


n 07+0-11 

u - y ' -0.18 

i-04±S:Si 
o.57±g : g 

n 25+ - 21 


0-57±8:l| (56) d 
0.62lg;S (56) 
0.58lg;£ (55) 

('.in (55) 
o.o8i8:8i (53) 


-1.40 
-1.31 
-1.29 
-1.54 
-2.09 


0.65±0.10 
0.79±0.33 
0.70±0.28 
0.41±0.07 
0.50±0.03 


0.94±0.14 
0.74±0.61 
0.94±0.65 
1.72±0.31 
1.04±0.18 



76+ 013 

u - ' °-0.23 


>1.053 


61 + 015 


(56) 


-1.29 


0.62±0.07 


0.57±0.10 


o.69±8:« 


0.954-1.053 


o.6ii8:iI 


(56) 


-1.31 


0.33±0.10 


0.99±0.18 


0.8518^ 


0.733-0.954 


0.6018:^ 


(55) 


-1.36 


0.41±0.02 


0.52±0.06 


61 +0 ' 32 


0.381-0.733 


o.i6l8:^ 


(56) 


-1.60 


0.47±0.01 


0.52±0.06 


^7+ - 31 
u - d ' -0.30 


<0.381 


o.o8i8:8§ 


(52) 


-2.07 


0.44±0.02 


0.23±0.04 



a We sort the spectra by the continuum slope and divide the objects into five groups, in which 
each group have the same number of sources (56). 

b The upper limits are also included. 

c The PAH luminosity is the sum of the 6.2, 7.7 and 11.3 //m PAH luminosities measured from 
the composite spectra. 

d The number in the parenthesis indicates the number of objects that we have measured the 
6.2 fim PAH EW or PAH luminosity. 
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5 10 15 20 25 5 10 15 20 25 

Wavelength (/xm) Wavelength (/u.m) 



Fig. 1. — Examples of PAHFIT decomposition of 5MUSES spectra (black squares) with 
strong, moderate, weak PAH emission and with silicate absorption. The best fit SED (green) 
is composed of thermal dust continua (red), PAH features (blue), stellar light (magenta) and 
emission lines (purple). 



-37- 



jrces 


60 
50 






J 


- 




5.8/xm 

rest — frame ; 


—> 
o 

00 


40 














'o 


30 














Number 


20 

10 





. _ 


i 









-0.4 -0.2 0.0 0.2 0.4 
log [(F(source)/F(tem plate)] 



w 



60 












8.0/xm 1 


o 
i_ 


50 












rest — frame - : 


o 


40 














in 














"o 


30 














CD 
_Q 


20 














E 
















D 


10 





1 








z 




1 

1 









-0.4 -0.2 0.0 0.2 0.4 
log [(F(source)/F(tem plate)] 



S 80 








24/xm 


u 








obs-frame 


o 60 

GO 










° 40 










<D 

E 20 
z 













-0.4 -0.2 0.0 0.2 0.4 
log [(F(source)/F(tem plate)] 




-0.4 -0.2 0.0 0.2 0.4 
log [F(source)/F(tem plate)] 



o 

00 



0) 
XI 

E 




-0.6 -0.4 -0.2 -0.0 0.2 0.4 0.6 
log [F(source)/F(tem plate)] 



Fig. 2. — The histogram of the flux ratios between the source and the best-fit template, for 
galaxies with MIPS 70 and 160 /zm detections (solid line) and for galaxies only with MIPS 
70 fim detections (dashed line). The dotted and dash-dotted lines indicate the medians of the 
flux ratios for sources with 70 and 160 /im detections and sources with only 70 fim detections. 
All five bands appear to peak around a log[(F(source)/F(template)] ratio of with rather 
narrow distributions. The la deviations in F (source) /F (template) for the 5.8, 8.0, 24 and 
70 160 /im bands are 0.07, 0.07, 0.03, 0.06 and 0.10 dex respectively for sources with 70 and 
160 fim detections and 0.07, 0.07, 0.03 and 0.07 dex for sources with only 70 fim detections. 
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Fig. 3. — SEDs of a sample of nine 5MUSES sources. The final SED (black) is composed of 
the IRS spectrum in the mid-IR and the best-fit template in the FIR. The observed data are 
shown as red circles. The three sources in the top panel are fit with 5 data points (IRAC 
5.8, 8.0 /mi and MIPS 24, 70 and 160 /im). The three sources in the middle panel are fit 
with 4 data points (IRAC 5.8, 8.0 /im and MIPS 24 and 70 /im). The three sources in the 
bottom panel are fit with only the IRS spectra. The blue line is the mid-IR SED of the 
best-fit template for the sources fit with the IRS spectra. 
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Fig. 4. — a) Top panel: The distribution of the total infrared luminosity of the 5MUSES 
sample. The dotted line represents the whole sample with known redshift. The blue, yellow 
and red solid lines represent the SB, composite and AGN sources in the sample. The SBs 
and AGN dominate the lower and higher end of the luminosity distribution respectively, b) 
Bottom panel: The distribution of the redshifts of the 5MUSES sample. The symbols are 
the same as in a). 
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Fig. 6. — The luminosity ratio of several bands over Lfr versus Lir. The blue circles, 
yellow crosses and red diamonds represent the SB, composite and AGN dominated sources 
in 5MUSES. The dotted, dash, and dash-dotted lines stand for the median ratios for the 
SB, composite and AGN sources, respectively. The PAH luminosities are derived from 
PAHFIT measurements. The 5.8, 14 and 24 fim luminosities are monochromatic luminosities 
calculated from the continua at these wavelengths. The IRAC 8 /im luminosities are derived 
by convolving the rest-frame IRS spectra with the filter response curve of the IRAC 8 /im 
band. The ratios and the associated uncertainties are also listed in Table [3j 
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Fig. 7. — a) Left panel: The median IRS spectra for SB, composite and AGN dominated 
sources from 5MUSES after normalizing at 5.8 /jm. The SEDs have been offset vertically. 
The shaded regions represent la uncertainties, b) Right panel: The median SEDs for SB, 
composite and AGN sources from mid-IR to FIR, normalized at 5.8 /im. 
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Fig. 8. — The histogram shows the distribution of the 6.2 /xm PAH EW for the 280 known-z 
sources from the 5MUSES sample. The solid line represents sources which have detection 
for the 6.2 /xm feature, while the dotted line includes upper limits. It is clear that both the 
solid line and dotted line show a dip in the PAH EW distribution at 0.2~0.3 fxm. See the 
text for detailed discussion on the bi-modality of the distribution. 
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Fig. 9. — a) Left panel: The continuum flux ratio of f7o Mm /f24/xm versus the 6.2 /zm PAH EW. 
b) Right panel: The continuum flux ratio of honm/hsum versus the 6.2 /zm PAH EWs. The 
solid line is a fit to all the data points while the dashed line is a fit excluding sources with 
6.2 /zm PAH EW upper limits. We bin the objects according to their continuum slopes and 
have equal numbers of objects in each bin. The dotted lines indicate the boundaries of those 
bins. 




Fig. 10. — a) Left panel: The typical mid-IR SEDs in each bin of different f 70^^24^111 ratios, 
b) Right panel: The typical SEDs in each bin of different hofim/hsttm ratios. All the SEDs 
have been normalized at 5.8 /im. The colors represent the average spectra derived in each 
ho^m/h^m (or f 3 0Mm/fi5^m) color bins listed in Table H 
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Fig. 11. — a) Left panel: The ratios of LiR/L 2 4 Mm versus L7o Aim /L 2 4 Aim - The sources are 
colored according to their fyo/f24 ratios. The dotted line is a 2nd-order polynomial fit to the 
data, b) Right panel: The ratios of L IR /L 15Mm versus L 30Mm /L 15Atm . The sources are colored 
according to their f3o/fis ratios. The dotted line is a linear fit to the data. 
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Fig. 12. — The PAH luminosity ratio of LpAH7.7/mi/LpAHii.3^m versus the 6.2 /mi PAH EW 
for 5MUSES. The AGN-dominated sources appear to have lower LpAH7.7/im/LpAmi.3/im ratios 
than the composite or SB dominated sources. The mean ratios are 3.45±0.55, 3.65±1.28 
and 2.26±0.89 respectively for SB, composite and AGN. 
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Fig. 13. — a) Top panel: A comparison of the distribution of PAH luminosity ratio of 
LpAH7.7pim/LpAHii.3/im for the IR classified SB-dominated sources from 5MUSES, SINGS 
and the optically classified SF-dominated sources in SSGSS. b) Bottom panel: Same as 
a), but for the PAH luminosity ratio of LpAH6.2/jm/LpAH7.7Aim- The SB and SF galaxies in 
5MUSES and SSGSS appear to have a similar distribution for both the LpAH7.7/mi/LpAHii.3/mi 
and LpAH6.2/mi/LpAH7.7/jm ratios, while the SINGS nuclear spectra appear to show higher 
LpAH7.7/jm/LpAHii.3/jm ratios for the SB-dominated galaxies. 
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Fig. 14.— The PAH band-to-band ratios of Lpahii 3/i.tti/Lpatt7.7/i.tti versus 

LpAH6.2/im/LpAH7.7/^m- The lines represent the expected ratios for neutral (upper line) 
and ionized (lower line) PAHs from model predictions. Note we have fewer objects in this 
Figure than in Figure [12], because we require the source to have S/N>3 for all three PAH 
features (6.2, 7.7 and 11.3 /xm) to be included. 
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Fig. 15. — The flux ratios of [NeIII]/[NeII] versus the PAH luminosity ratios of 
LpAH7.7/im/LpAHii.3/im- The blue, yellow and red open triangles represent upper/lower limits 
for the SB, composite and AGN in 5MUSES and the directions the triangles face are con- 
sistent with the directions of the limits. The dotted line is the median LpAH7.7/xm/LpAHii.3/mi 
ratio for the SB-dominated sources on this plot. 



